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Abstract 

The Henon-Heiles system in the general form has been considered. In a 
nonintegrable case with the help of the Painleve test new solutions have been 
found as formal Laurent or Puiseux series, depending on three parameters. One 
of parameters determines a location of the singularity point, other parameters 
determine coefficients of series. It has been proved, that if absolute values of 
these two parameters are less or equal to unit, then obtained series converge 
in some ring. For some values of these parameters the obtained Laurent series 
coincide with the Laurent series of the known exact solutions. 



1 THE PAINLEVE PROPERTY AND INTEGRABI- 
LITY 

A Hamiltonian system in a 2s-dimensional phase space is called completely integrable (Li- 
ouville integrable) if it possesses s independent integrals which commute with respect to the 
associated Poisson bracket. When this is the case, the equations of motion are (in principal, 
at least) separable and solutions can be obtained by the method of quadratures. 

When we study some mechanical or field theory problem, we imply that time and space 
coordinates are real, whereas the integrability of motion equations is connected with the 
behavior of their solutions as functions of complex time and (in the case of the field theory) 
complex spatial coordinates. 

S.V. Kovalevskaya was the first, who proposed Jl| to consider time as a complex variable 
and to demand that solutions of the motion equations have to be single- valued, meromorphic 
functions on the whole complex (time) plane. This idea gave a remarkable result: S.V. Ko- 
valevskaya discovered a new integrable case (nowadays known as the Kovalevskaya 's case) 
for the motion of a heavy rigid body about a fixed point |l| (see also [fj, ||). The work of 
S.V. Kovalevskaya has shown the importance of application of the analytical theory of differ- 
ential equations to physical problems. The essential stage of development of this theory was 
a classification of ordinary differential equations (ODE's) in order of types of singularities of 
their solutions. This classification has been made by P. Painleve. 

Let us formulate the Painleve property for ODE's. Solutions of a system of ODE's 
are regarded as analytic functions, may be with isolated singular points ||, |5|. A singular 
point of a solution is said critical (as opposed to noncritical) if the solution is multivalued 
(single- valued) in its neighborhood and movable if its location depends on initial conditions^. 

Definition. A system of ODE's has the Painleve property if its general solution 
has no movable critical singular point |6|. 0]. 

An arbitrary solution of such system is single-valued in the neighborhood of its singular 
point t and can be expressed as a Laurent series with a finite number of terms with negative 
powers of t— 1 . If a system has not the Painleve property, but, after some change of variables, 
the obtained system possesses this property, then the initial system is said to have the weak 
Painleve property. 

Investigations of many dynamical systems, Hamiltonian [8-10] or dissipative (for example, 
the Lorenz systems [10-13]), show, that a system is completely integrable only for such 
values of parameters, at which it has the Painleve property (or the weak Painleve property). 
Arguments, which clarify the connection between the Painleve analysis and the existence of 
motion integrals, are presented in [[TJ], [OJ. If the system misses the Painleve property (has 



complex or irrational "resonances"), then the system cannot be "algebraically integrable" [16| 



(see also |T7| and references therein). At the same time the integrability of an arbitrary 
system with the Painleve property has yet to be proved. There is not an algorithm for 
construction of the additional integral by the Painleve analysis. It is easy to give an example 



Solutions of a system with a time-independent Hamiltonian can have only movable singularities. 
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of an integrable system without the Painleve property [18|: H = ^p 2 + f(x), where f(x) is a 
polynomial which power is not lower than five. The given system is trivially integrable, but 
its general solution is not a meromorphic function. 

The Painleve test is any algorithm designed to determine necessary conditions for a 
differential equation to have the Painleve property. The original algorithm, developed by 
P. Painleve and used by him to find all the second order ODE's with Painleve property [0, 
is known as the a-method. The method of S.V. Kovalevskaya is not as general as the 
CK-method, but much more simple^. 

In 1980, motivated by the work of S.V. Kovalevskaya [|1], M.J. Ablowitz, A. Ramani and 
H. Segur |2(J developed a new algorithm of the Painleve test for ODE's. The remarkable 
property of this test is that it can be checked in a finite number of steps. They also were 
the first to point out the connection between the nonlinear partial differential equations 
(PDE's), which are soluble by the inverse scattering transform method, and ODE's with 
the Painleve property. Subsequently the Painleve property for PDE was defined and the 
corresponding Painleve test (the WTC procedure) 



(see also 19 



21 



was constructed 

23-29]). With the help of this test it has been found, that all PDE's, which are solvable 
by the inverse scattering transforms, have the Painleve property, may be, after some change 
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of variables. For many integrable PDE's, for example, the Korteweg-de-Vries 
sine-Gordon [2j| equations, the Backlund transformations and the Lax representations result 
from the WTC procedure [22, 21, |29|. Also, special solutions for certain nonintegrable PDE's 
were constructed using this algorithm [|3(], [31|| . 

The algorithm for finding special solutions for ODE system in the form of a finite ex- 
pansion in powers of unknown function ip(t — to) has been constructed in |32, 33]. The 
function (p(t — to) and coefficients have to satisfy some system of ODE, often more simple 
than an initial one. This method has been used |34| to construct exact solutions for certain 
nonintegrable systems of ODE's. 

The aim of this paper is to find new special solutions for the generalized Henon-Heiles 
system using the Painleve test. In distinction to [34] we obtain solutions as formal Laurent 
or Puiseux series and find domains of their convergence. 



2 THE HENON-HEILES HAMILTONIAN 

Let us consider the motion of a star in an axial-symmetric and time-independent potential. 
The motion equations admit two well-known integrals (energy and angular momentum) and 
would be solved by the method of quadratures if the third integral of motion is known. Due 
to the symmetry of the potential the considered system is equivalent to two-dimensional one. 
However, for many polynomial potentials the obtained system has not the second integral 
polynomial function. 

In the 1960s, asymptotic methods |fJ3|. have been developed to show either existence 
or absence of the third integral for some polynomial potentials. To answer the question 
about the existence of the third integral Henon and Heiles |37| considered the behavior 
of numerically integrated trajectories. Emphasizing that their choice of potential does not 
proceed from experimental data, Henon and Heiles have proposed the following Hamiltonian: 

2 Different variants of the Painleve test are compared in pi R. Conte paper] 
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H = - (x 2 + y 2 + x 2 + y 2 ) + x 2 y - ^y 3 , (1) 

because: on the one hand, it is analytically simple; this makes the numerical computations 
of trajectories easy; on the other hand, it is sufficiently complicated to give trajectories 
which are far from trivial. Indeed, for low energies the Henon-Heiles system appears to be 
integrable, in so much as trajectories (numerically integrated) always lay on well-defined two- 
dimensional surfaces. On the other hand, for high energies many of these integral surfaces 
are destroyed, it points on absence of the third integral. 



Subsequent numerical investigations fl38|, |39| show, that in the complex t-plane singular 




points of solutions of the motion equations group in self-similar spirals. It turns out extremely 
complicated distributions of singularities, forming a boundary, across which the solutions can 
not be analytically continued. 

The generalized Henon-Heiles system is described by the Hamiltonian: 

H = \{x 2 t + y 2 + A* 2 + y 2 ) + x 2 y-j y 3 (1') 
and the corresponding system of the motion equations: 

= —Ax — 2xy, 

(2) 

= -y-x 2 + Cy 2 , 

where x u = and y tt = ^J[, A and C are numerical parameters. 

Due to the Painleve analysis the following integrable cases of (2) have been found: 

(i) C=-l, A = l, 

(ii) C = —6, A is an arbitrary number, 

(iii) C=-16, A = i. 

In contradiction to the case (i) the cases (ii) and (iii) are nontrivial, so the integrability 
of these cases had to be proved additionally. In the 1980 's the required second integrals 
were constructed [40-43]. For integrable cases of the Henon-Heiles system the Backlund 
transformations pi R31 and the Lax representations P5], |2"E|, [4"5| have been found. 



The three integrable cases of the Henon-Heiles system correspond precisely to the sta- 
tionary flows of the only three integrable cases of fifth- order polynomial nonlinear evolution 
equations of scale weight 7 (respectively the Sawada-Kotega, the fifth-order Korteweg-de 



Vries and the Kaup-Kupershmidt equations) |53], [44], [45 1. 

The Henon-Heiles system is a model widely used in physics, in particular, in gravita- 
tion [46-48] and plasma theory |49]. The models, described by the Hamiltonian (!') with 



some additional nonpolynomial terms, are actively studied HR 50-52] as well. 



3 NONINTEGRABLE CASES 



The general solutions of the Henon-Heiles system are known only in integrable cases |52j| , in 
other cases search of new (exact or asymptotic) solutions is an actual problem. 
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The procedure for transformation the Hamiltonian to a normal form and for construction 
the second independent integral in the form of formal power series in the phase variables 
x, x t , y and y t (Gustavson integral) has been realized for the Henon-Heiles system both in 
the original (A = 1, C — 1) |3B|] (see also pBfl ) and in the general forms |>|, Using the 
Bruno algorithm [|56], |57| V.F. Edneral has constructed the Poincare-Dulac normal form 
and found [59| (provided that all phase variables are small) local families of periodic 
solutions. Recently it has been found that a local series around the singularities in the 
complex (time) plane can be transformed to some local series around the singularities at the 
fixed points in phase space and analyzed via normal forms theory |K| |H[ . 

The Henon-Heiles system as a system of two second order ODE's is equivalent to the 
fourth order equationP]: 

y tttt = (2C - 8)y tt y - (4A + l)y tt + 2{C + l)y 2 t + frac20C3y 3 + (4CA - 6)y 2 - \y - AH, (3) 

where H is the energy of the system. 

To find a special solution of the given equation one can assume that y satisfies some more 
simple equation. For example, the well-known solutions in terms of the Weierstrass elliptic 
functions |62], 



x| satisfy the following first-order differential equation: 
y 2 = Ay 3 + By 2 + Cy + V, 



where 



or 



A 



B 



C 



V 



A = -C, B = -1, C = and V = 2H 



4 
3 : 



(4) 



(4a) 



1 - (C + 2)A 



C + l 
3C 2 X 2 - 



3C 2 A + 8CA 2 - 7CX - C + 4A 2 - 2A - 2 



3C 3 + IOC 2 + 11C + 4 
2AC 4 H + 104C 3 H - 9C 3 X 3 + 6C 3 X 2 + 3C 3 X 



(46) 



+ 



+ 



4(3C 5 + 22C 4 + 60C 3 + 78C 2 + 49C + 12) 

168C 2 H - 30C 2 X 3 + 13C 2 X 2 + 16C 2 X + C 2 
4(3C 5 + 22C 4 + 60C 3 + 78C 2 + 49C + 12) 

120CH - 28CX 3 + 2ACX + AC + 32H - 8X 3 



+ 



+ 



4A 2 + 8A + 4 



4(3C 5 + 22C 4 + 60C 3 + 78C 2 + 49C + 12) 

T> is proportional to energy H (arbitrary parameter), therefore, solutions (4a) and (46) are 
two-parameter ones. 



3 For given y(t) the function x 2 (t) is a solution of a linear equation. System (2) is invariant to exchange 
x to —x. 
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E.I. Timoshkova [64] generalized equation (4): 



Ay 3 + By 2 + Cy + V + Gy 5/2 + £y 3 ' 2 (5) 



and found new one-parameter sets of solutions of the Henon-Heiles system in nonintegrable 
cases (C = — | or C = — A is an arbitrary number). These solutions (i.e. solutions 
with Q ^ or £ ^ 0) are derived only at V = 0, therefore, substitution y = g 2 gives: 

g 2 = ^(Ag 4 + gg 3 + Bg 2 + £g + Cy (6) 

The general solution of (6) has one arbitrary parameter and can be expressed in elliptic 
functions. 

In this paper I analyze system (2) at C = — and A = | (the Solution 2.2 of the 
paper [J64]] ). In this case equation (5) is: 

and, depending on a choice of a sign before the last term, we obtain either (in case of sign 
+): 

K 25(l- sin (^)) 

and x 2 = — ^ 4; (8.1) 



3(l-3sin(^)) 9 (1 -3sin(^a)) > 

or (in case of sign — ): 

R 25fl + sin (^a) 

and x 2 = — ^ (8.2) 



3(l + 3sin(^)) 9 (1 + 3sin(^))' 

4 RESULTS OF THE PAINLEVE TEST FOR THE 
HENON-HEILES SYSTEM 

The Ablowitz-Ramani-Segur algorithm of the Painleve test appears very useful to find 
asymptotic solutions as a formal Laurent series. 

We assume that the behavior of solutions in a sufficiently small neighborhood of the 
singularity is algebraic, it means that x and y tend to infinity as some powers of t — to: 

x = a a (t — t ) a and y = bp(t — to)' 3 ; (9) 

where a, /3, a a and bp are some constants. We assume that real parts of a and (3 are less 
then zero, and, of course, a a 7^ and bp 7^ 0. 

If a and (3 are integer numbers, then substituting 



* max 



N„ 

x = a a (t - t ) a + °i+a(* - *o) i+Q , (10.1) 

3=1 
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Nmax 



3+P 



'10.2) 



3=1 



one can transform the ODE system into a set of linear algebraic systems in coefficients 
and In the general case one can obtain the exact solutions (in the form of formal Laurent 
series) only if one solves infinity number of systems (N max = oo). On the other hand, if one 
solves a finite number of systems one obtains asymptotic solutions. With the help of some 
computer algebra system, for example, the system REDUCE [^5[ |66|, these systems can be 
solved step by step and asymptotic solutions can be automatically found with any accuracy. 
But previously one has to determine values of constants a, f3, a a and bp and to analyze 
systems with zero determinants. Such systems correspond to new arbitrary constants or 
have no solutions. Powers at which new arbitrary constants enter are called resonances. The 
Painleve test gives all information about possible dominant behaviors and resonances (see, 
for example, [TIJ). Moreover, the results of the Painleve analysis point out cases, in which 
it is useful to include into expansion terms with fractional powers of t — t . 

For the generalized Henon-Heiles system there exist two possible dominant behaviors and 
resonance structures |I| pSfl : 



Case 1: 


Case 2: {(3 < 9fce(a) < 0) 


a = -2, 
P = -2, 

a a = ±3^2 + C, 
bp = -3, 


\±Ji-m/c 

ry = i 

2 ' 

(3 = -2/ 

a>a = C\ (an arbitrary number), 

b - ^ 

°p — c" 


r _ i fi 5 . y/1-up.+C) 

' ' 2 2 


r = -l, 0, 6, Tv /l-48/C. 



In the Table the values of r denote resonances: r = — 1 corresponds to arbitrary parameter 
t Q ; r = (in the Case 2) corresponds to arbitrary parameter c 1 . Other values of r determine 
powers of t, to be exact, t a+r for x and t^ +r for y, at which new arbitrary parameters enter 
(as solutions of systems with zero determinants). 

For integrability of system (2) all values of a and r have to be integer (or rational) and 
all systems with zero determinants have to have solutions at all values of included in them 
free parameters. It is possible only in the cases (i) — (iii). 

At C = —2 (in the Case 1) a a = 0. It is the consequence of the fact that, contrary to 
our assumption, the behaviour of the solution in the neighborhood of a singular point is not 
algebraic, because its dominant term includes logarithm. 

Those values of C, at which a and r are integer (or rational) numbers either only in the 
Case 1 or only in the Case 2, are of interest for search of special solutions. 
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5 NEW SOLUTIONS 



5.1 Finding of solutions in the form of formal Laurent series 

Let us consider the Henon-Heiles system with C = — ^f. In the Case 1 some values of r 
are not rational, so it is a nonintegrable system. To find special asymptotic solutions let us 
consider the Case 2. In this case a — — § and r = —1, 0, 4, 6, hence, in the neighborhood 
of the singular point to we have to seek x in such form that x 2 can be expand into Laurent 
series, beginning from (t — to) -3 - Let t = 0, substituting 



x 



\ 15 

= Vt ( c x t~ 2 + a i tj I and y = — 1~ 2 + M j 

J I / 1 j=-l 



in (2), we obtain the following sequence of linear system in a k and b k : 



k-i 



(k 2 - A s ja k + 2c ± b k = -Aa fc _ 2 - 2 ^ ajb k -j- 2 

j=-i 

k-l ^ fe-1 

({k - l)k - 12^)b k = - b k - 2 ~ Y a 3 a k-j-i ~ y Y b J b k-j~2- 



(11) 



5 

j=-2 J=-l 



If /c = 2 or A; = 4, then the determinant of (11) is equal to zero. To determine a 2 and b 2 
we have the following system: 



' q (557056c? + (15552000A - 4860000)ct + 8640000006 2 + 
+ 108000000A 2 - 67500000A + 10546875) = 0, 

818176c? + (15660000A - 4893750^- 
[ - 810000000&2 - 6328125 = 0. 



(12) 



As one can see this system does not include terms, which are proportional to a 2 , hence, 
a 2 is an arbitrary parameter (a constant of integration). 

We discard the solution with c x — and obtain the system in c x = c\ and b 2 with the 
following solutions: 



Ci = 



1125(4 V /35(2048A 2 - 1280A + 387) - 1680A + 525) 

167552 ' 
(10944A - 3420) v/35(2048A 2 - 1280A + 387) - 4403456A 2 + 2752160A - 789065 
2 117956608 

or 

1125(-4^/35(2048A 2 - 1280A + 387) - 1680A + 525) 
Cl ~ 167552 ' 

(10944A - 3420) v/35(2048A 2 - 1280A + 387) - 4403456A 2 + 2752160A - 789065 
2 ~ 117956608 ' 



We obtain new constant of integration a 2 , but we must fix c 1; so number of constants 
of integration is equal to 2. It is easy to verify that 6 4 is an arbitrary parameter, because 
the corresponding system is equivalent to one linear equation. So, using Painleve test, we 
obtain an asymptotic solution which depends on three parameters, namely t , a 2 and 64. 

Now asymptotic solutions can be obtained with arbitrary accuracy For given A one has 
to choose c l as one of the roots of system (12). After this the coefficients aj and bj can be 
found automatically with the help of some computer algebra system. 

For example, if A = |, then (12) has the following solutions: 

_ 625 1819 1 f 8125 8700683 1 

L ~ 128' 2 " " 663552 J ' \ Cl " ~ 23936' 2 " " 1364926464 J ' 

Taking into account, that system (2) is invariant to change x to —x, we obtain four types 
of formal solutions: 

r \h^2 2 25 1 5^2 5275 2 

x = Vt < — — t" 2 + — -^r 1 — + —t + a 2 t 2 

\ 4 96^2 16 663552^2 

15 2 5^2 , 205 115V2 1819 2 no in 

y = r 2 + -^— r 1 + — —t t 2 + 13.1 

y 8 32 2304 13824 663552 

/ 741719^ 5^2 . ., . 

+ h — —a 2 t 3 + b 4 r + ...: 

^1528823808 12 ! 



x 



r \^i</2 2 25i , 5i^2 5275i , 2 

Vt < — — r 2 -r 1 -t 1 + a 2 t 2 + 

1 4 96^2 9216 663552^2 



15 2 5v/2 , 205 115v^ 1819 2 /„m 

w = — t" 1 — t 1 2 - (13.2) 

y 8 32 2304 13824 663552 

/ 741719^ 5iv^2 , . 
" ^1528823808 + ' ^ +U + -> 



.7; 



1 4 V 374 17952 V 374 



38645^ J 13~ 7028575i v / 243l J 13~ 2 
~ 574464 V ~ 374 ~ 23203749888 V ~ 374* + + ' ' ' 8 ' ' 

15 _ 2 hiV^Ml _! 69335 37745iv / 4862 _ 8700683 
^ ~ " T f 5984 ~ 430848 ~~ 483411456 1364926464 



\ 3332429743915008 12 V 374 



/ 1148020763iv / 13v / 374 5^2 J 13" , ,, , , 



13.3) 
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r , 5i\/2 A 13" _ 2 25 V / 243T A 13" 



38645z^2 A 13" 7028575 v / 2431 A 13~ 2 
574464 V ~ 374 ~ 23203749888 V _ 374* + ° 2 

13.4) 

15 _ 2 5«V4862 _j 69335 37745^4862 8700683 2 
~ ~ ~8~ 5984 * ~ 430848 + 483411456 * ~ 1364926464* 

/ 1148020763/0^374 , 5i^2 4 / 13" , ,, . , , 



+ -r^^2\ ~ 7737 ^ + &# + •••■ 



V 3332429743915008 12 V 374 
It is easy to verify that if 



21497^2 , , 858455 

0,2 = — - and 



42467328 12039487488 
then series (13.1) are the Laurent series of (8.1). Also if 

21497^^2 , , 858455 

a>2 = and 



42467328 12039487488 
then series (13.2) are the Laurent series of (8.2). 

5.2 Convergence of the obtained series 

When an asymptotic series is obtained the question about its convergence arises. It is known 
that a domain of Laurent series convergence is a ring. Let us find conditions, at which the 
obtained series converge in the following ring: < \t\ ^ 1 — e, where e is any positive number. 

The sum of a geometrical progression S = t n = is finite if \t\ ^ 1 — e, hence, our 

n=0 

series converge in the above-mentioned ring, if 3N such that Vn > N \a n \ ^ 1 and \b n \ ^ 1. 
Let \a n \ ^ 1 and \b n \ ^ 1 for all — 1 ^ n < k, then from (11), we obtain: 

. 2/c + 2+ |A| +2| Cl | . 21(k + 2) 

|afc| < \k^A\ ' W < 5(P- fe -12) - (M) 

It is easy to see that there exists such N, that, if \a n \ ^ 1 and \h n \ ^ 1 for — 1 ^ n ^ N, 
then \a n \ ^ 1 and \b n \ 1 for — 1 ^ n < 00. N is a maximum from 8 and l + -\/|A|+2|c 1 | + 7. 

For example, if A = |, then for any possible value of c l3 we obtain N = 8. It is easy to 
verify that if |ci 2 1 ^ 1 and 1 1> 4 1 ^ 1, then \a n \ ^ 1 and \b n \ ^ 1 for — 1 ^ n ^ 8, and, hence, 
for an arbitrary n. Thus our Laurent series converge in the ring < \t\ ^ 1 — e. Numerical 
analysis shows |H| that these series can also converge at absolute values of parameters more 
than unit. For other values of A the convergence can be considered analogously. 
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6 GENERALIZATION OF SOLUTIONS IN TERMS 
OF THE WEIERSTRASS ELLIPTIC FUNCTIONS 

Let us consider solutions (4). The values of parameters (4a) correspond to x(t) = 0. Solu- 
tions (46) correspond to the Case 1 (see Table) and x(t) can be expressed in terms of the 
Weierstrass elliptic functions. For some values of C these two-parameter solutions can be 
generalized. For example, if C — — |, then some resonances are half-integer. Substituting 



oo 
fc=-4 



3V7+-2 , _J_ , _±_ n -3/2 , V7 + 2 _ 15y/7 n f 7/2 _ 
2>/2 C + Sv^ + 1 + 160^ 224v^ 1 



in (2), we obtain that a n and b n have to satisfy the following system: 

n(n-2) -24_ ~ ^ _ ~ 
a n + 2a_40 n — — Aa„_4 — 2 akb n -k-4, 

k=-i 

o~ ~ n{n-2)-TI~ ~ ^ _ „ 9 ^ ~ ~ 

2a_4d„ H n — — n _ 4 — a k0. n -k-1 ~ g ^ ^kbn-k-A, 

k=-3 k=-3 

where a_4 = ± and o_ 4 = —3. 

At any A we obtain three-parameter solutions as formal Puiseux series. For example, if 
A = 1 (and a_ 4 = |^|) the solution is the following (t = 0): 

x = 

_ _vfL D +4 _ 467^ n 2 f 5 1157^7 n ,11/2 , y 7 ? -6 , 

2v^ 2 8624^ 1 + 430080^ 1 115200>/2 + ' * * 

y = -3r 2 -| + A^ 2 -^ 2 -i^ 2 + 

, n 79 n2 f 5 I 1157 n .11/2 1 y.6 i 

+ U 2 t - eTe^lI + 245760-^ 57600 t + " " " " 

Using numerical calculations it is easy to show that if < 1 and \D 2 \ < 1 then \hk\ < 1 

and \b~k\ < 1) for — 3 ^ k ^ 50, except only a 3 = ^J^D X . It is sufficient to prove that \a k \ < 1 

and \bk\ < 1 for all k > 50 and, hence, our series converge in the ring < \t\ ^ 1 — e. If 
Dx — then y satisfies (4) with 

A=--, B=-l, C = and V=—H, (46') 
3 15 

and solution can be presented in terms of the Weierstrass elliptic functions. 



(15) 



7 CONCLUSION 

Using the Painleve analysis one can not only find integrable cases of dynamical systems, but 
also construct special solutions in nonintegrable cases. 



11 



We have found the special solutions of the Henon-Heiles system with C = - y as for- 
mal Laurent series, depending on three parameters. For some values of two parameters the 
obtained solutions coincide with the known exact solutions. At C — — | two-parameter solu- 
tions in terms of the Weierstrass elliptic functions have been generalized to three-parameter 
ones. New solutions found as formal Puiseux series. For some values of A the analysis of 
convergence of the obtained series has been made and it has been proved, that they have 
nonzero domain of convergence. Similar analysis can been made for any value of A. 

With the help of the Painleve test particular asymptotic solutions as Laurent or Puiseux 
series can be found for the Henon-Heiles system with some other values of C and A. Just 
at these values of parameters the probability of finding of new exact solutions similar to the 
solutions found in |63| is great. 

The author is grateful to R. I. Bogdanov and V. F. Edneral for valuable discussions 
and E. I. Timoshkova for comprehensive commentary of ^M. This work has been supported 
by the Russian Foundation for Basic Research under grants A/° 00-15-96560 and 00-15-96577 
and by the program " Universities of Russia" . 
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